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Abstract 

£> ' Using the newest data for pp scattering at the CERN Large Hadron Collider (LHC) 

combined with the Glauber model, we make hydrodynamic predictions for the soft 
pg \ hadronic observables planned to be measured in the forthcoming Pb+Pb collisions 

(N ; at v /s^i = 2.76 TeV. 

PACS: 25.75.-q, 25.75.Dw, 21.65.Qr 

1. The behavior of matter created in relativistic heavy-ion collisions at the 
BNL Relativistic Heavy Ion Collider (RHIC) is very well described by the 
perfect-fluid hydrodynamics [1.2.3.4.5.6.7.8.9]. Having in mind these successes, 
it is appealing to use perfect-fluid hydrodynamics to make predictions for the 
Pb+Pb collisions at the LHC. Several such studies have been performed so 
far and they are summarized in Ref. [TO] . 

Nevertheless, in the last two years two important aspects concerning the 
heavy-ion physics at the LHC have changed: Firstly, the forthcoming heavy-ion 
experiments will have a lower energy than that planned previously, a /snn = 2.76 
TeV instead of ^/snn = 5.5 TeV. Secondly, the data for pp scattering at the 
LHC have just been analyzed. These new facts suggest that a revised predic- 
tion for ^/snn = 2.76 TeV based on the newest pp results is advisable. 

In the paper we follow this idea and use the newest data for pp scattering 
at the LHC to make hydrodynamic predictions for soft hadronic observables 
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measured in Pb+Pb collisions at ^/snn = 2.76 TeV. The pp data are rescaled 
by the number of sources in heavy-ion collisions (estimated from the Glauber 
model) in order to determine the final multiplicity of charged hadrons in nu- 
clear collisions. The relative increase of the multiplicity of particles produced 
at RHIC and LHC energies requires a similar increase of the initial entropies. 
This, in turn, allows us to make estimates of the initial central temperature 
expected at the LHC. 

Our hydrodynamic approach is defined in Ref. [8], where we have described 
quite successfully the RHIC data and made predictions for Pb+Pb collisions 
at ^/snn = 5.5 TeV. It is a 2+1 dimensional (boost-invariant) hydrodynamic 
model with a realistic equation of state based on the lattice results and the 
hadron-gas model. For simplicity, a single-freeze-out scenario is adopted in the 
Monte-Carlo version [UJ - all known hadrons are generated on the freeze-out 
hypersurface determined by the condition of constant freeze-out temperature, 
then, unstable resonances decay with the branching ratios given by the PDG. 
The hadronic interactions in the final state are neglected. 

2. In order to specify the initial entropy density for our hydrodynamic cal- 
culations we determine first the expected multiplicity in nuclear collisions. In 
Ref. [12] the data on dN pp /dr/\ v=0 from the first proton-proton collisions at 
the LHC [T3fl4||15|ll6|ll7] have been fitted with the parameterization 
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where E is the CMS energy of the collision. For E = 2.76 TeV this fit gives 
dN pp /drj\ r]= o = 4.68. We shall, however, rather base our estimate on the closest 
data point which has been measured by ALICE [Qj and the CMS Collabora- 
tion [T7] at E — 2.36 TeV. The weighted average of the two results with weights 
given by the inverse of statistical uncertainties gives dN pp (2. 36) /drj \ v=0 = 4.65. 
This result is extrapolated with the power-law prescription flTJ to 
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4.65 = 4.82. 



This value will be used in our further calculations. 

The expected multiplicity in Pb+Pb collisions is obtained by the multiplica- 
tion of the pp result by the number of sources that are identified with the 
wounded (participant) nucleons and binary collisions 
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Here the number of participant nucleons N w , binary collisions N^ in , and the 
number of sources N ST depend on centrality and are obtained from the GLIS- 
SANDO code [T8] that is a Monte-Carlo implementation of the Glauber model. 
The parameter a in (|3j) takes into account the contribution of minijets to the 
fireball density. The jet contribution is expected to rise with the energy of 
the collisions. The centrality studies performed by PHOBOS on Au+Au col- 
lisions gave a = 0.12 at ^J]^ = 19.6 GeV and a = 0.145 at 200 GeV p3]. 
In the present analysis, assuming the linear growth of a with log(s), we take 
a = 0.16. 



The centrality dependence of the charged particle pseudorapidity density fol- 
lows from the centrality dependence of N SI (b) that is obtained from GLIS- 
SANDO with NN cross section of 63 mb. It is presented in Fig. [IJ 

3. For our hydrodynamic simulations we need the initial profiles of the entropy 
density. Following the standard approach used in hydrodynamic models we 
assume that the initial entropy density is proportional to the source density 
in the transverse plane 



s(b,x ± ) = s(T t ) 
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Fig. 1. Charged particle multiplicity density in pseudorapidity for Au+Au collisions 
at y^SNN = 200 GeV (solid line) together with the PHOBOS data [20], and predic- 
tions for Pb+Pb collisions at 2.76 TeV using a = 0.16. The model calculations are 
performed with GLISSANDO [E] and use Eqs. © and @. 



where s(T») is the central entropy density at zero impact parameter. The 
integral of the source density defines iV sr , 

N sr (b)= fd 2 x ±PsT (b,x ± ). (6) 



Unfortunately, GLISSANDO does not provide the source densities which are 
sufficiently smooth and could be directly used as initial conditions for our 
simulations. We thus calculate them in the optical Glauber approach in which 
we fix the impact parameter b in such a way that we reproduce the number 
of sources obtained in GLISSANDO, in agreement with Eq. (J6]). 

The absolute scale of the entropy density is estimated from the ratio of the 
multiplicity densities. We assume that the ratios of total entropies in the 
transverse plane (at midrapidity and at the same initial proper time) should 
be equal to the ratios of the multiplicities. By considering the ratios of the 
entropies for the LHC (at ^s^ = 2.76 TeV) and RHIC (at ^s^ = 200 GeV) 
we find 
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Expressing the multiplicity density at the LHC with the help of Eq. ([3]) we 
obtain the ratio of the entropy densities 

s(T PbPb ) /C P6 (0,0) 4.82 N AuAu (b) 
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The calculations of the source densities give 



p PbPb (0, 0) = 6.76 mi^, p AuAu (0, 0) = 4.71 far 2 . (9) 

For central Au+Au collisions at ^snn = 200 GeV (centrality class 0-5%) we 
find N AuAu = 308 and the PHOBOS multiplicity density is 667 (2D]. The ratio 
of entropies in Eq. ((7|) is then 3.2, and the ratio in Eq. (JS)) acquires practically 
the same value |_j. The transverse size of the fireball S± = tc J (x 2 ) (y 2 ) = 
22.9 fm 2 for the most central Au+Au collisions at RHIC is almost the same 
as for central Pb+Pb collisions at the LHC: Si = 23.2 fm 2 . 



1 For larger values of b, using the Glauber estimates for N A u and the experi- 
mental values for dN /dr] L=o> we fi n d deviations from 3.2 by about 10%. Such 
differences indicate that the scalings (TTJ and J5]) are valid with a similar accuracy. 
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Fig. 2. Predictions for transverse-momentum spectra of pions, kaons, and protons 
for the two centralities: c=0-5% (upper part) and c=20-30% (lower part). The 
values of the input parameters are discussed in the text. 

The entropy density is connected to the temperature through the equation of 
state constructed in [21J. The multiplicity in Au+Au collisions at ^snn = 200 
GeV is reproduced if the initial temperature at Tj = 0.25 fm is 520 MeV. From 
this we derive the initial temperature for ^/snn = 2.76 TeV to be 735 MeV 
at the initial proper time r, = 0.25 fm (the one used in the present paper), or 
480 MeV if n = 1 fm. 

4. In Figs. [2H31 we present our results for: the transverse-momentum spec- 
tra of pions, kaons, and protons, the elliptic flow of charged particles, and 
the HBT radii of pions. The calculations are done for two centrality classes 
(c = 0-5% and c = 20-30%). The hydrodynamic calculations use the follow- 
ing parameters: the starting time of the hydrodynamic evolution t% = 0.25 fm, 
the initial central temperature (in the most central collisions) Tj = 735 MeV, 
and the freeze-out temperature Tf = 150 MeV. From the measured p/p ratios 
in proton-proton collisions at the LHC [22] we estimate that p/p = 0.97 in 
nuclear collisions at y / J/v r /v = 2.76 TeV, which corresponds to a small value of 
the baryon chemical potential \xb — 2.2 MeV, if no additional baryon stop- 
ping occures in Pb+Pb collisions. In our calculations all chemical potentials 
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Fig. 3. The elliptic flow coefficient v% of charged particles shown as a function of the 
transverse momentum. The parameters the same as in the lower part of Fig. [2j 

in THERMINTOR have been set equal to zero. 

The spectra shown in Fig. [2] yield the following mean values of the transverse 
momentum: (pf ons ) = 656 MeV, ( v ^™) = 1028 MeV, and (pP rotons ) = 1460 
MeV for the centrality class c = 0-5% , and (p P i° ns ) = 637 MeV, ( p f°™) = 997 
MeV, and ( p P rotons ) = 1406 MeV for the centrality class c = 20-30%. The 
mean transverse momentum increases from the top RHIC energies to the LHC 
energies by about 35 — 40%. The expected pseudorapidity (rapidity) density 
of charged particles is dN/dr] = 2161 (dN/dy = 2457) for c = 0-5%, and 905 
(1036) for c = 20-30%. More specifically we predict 
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for 7r + , K + and protons respectively, for central collisions. The rise of the 
entropy from RHIC to LHC and the similarity of the transverse sizes of the 
fireballs leads to an increase of the upper range of the scaling density [23] from 
the value 
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at RHIC energies, to 
1 dN 



S± dy 



105 fm" 



(12) 



for the central Pb+Pb collisions at ^s^n = 2.76 GeV. 
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Fig. 4. The pionic HBT radii shown as functions of the transverse momentum of 
the pairs for the central (squares) and non-central (triangles) collisions. 

In Fig. [3] the elliptic-flow coefficient t>2 of charged particles is shown as a 
function of transverse-momentum. Compared to the analogous RHIC plots, 
the dependence of vi on p± exhibits saturation with increasing beam energy, 
as observed by Kestin and Heinz [21] . However, the substantial increase of the 
initial entropy density (matching the expected rapid growth of the multiplicity 
with energy) and the hard equation of state used in the calculations lead to 
a significant increase of (p±) and the integrated elliptic flow when going from 
from RHIC to the LHC. We find t>2 = 0.088 for charged particles in the 



centrality class c = 20 — 30%. The elliptic flow is known to be very sensitive to 
shear viscosity effects [25] that reduce its value. However, at the first stage of 
the expansion at the LHC, where stress corrections are the most important, we 
deal with high temperature and the relative importance of viscosity corrections 
is reduced compared to RHIC energies. 

FigureHJshows our predictions for the HBT radii. The radii are calculated with 
the methods defined in Ref . [26] . Here again, a small increase of the radii, as 
compared to RHIC, is observed. 

5. Summary: In this paper we make predictions for soft hadronic observables 
that will be soon measured in Pb+Pb collisions at a/snn = 2.76 TeV. Our 
results are based on the hydrodynamic model, which was successfully used to 
describe the RHIC results. Compared to the RHIC data, the values of various 
physical parameters show smooth changes with the beam energy, however, 
our predictions indicate larger changes of the multiplicity than those expected 
earlier. The main reason for this fact is the large multiplicity observed in pp 
collisions at the LHC, which is used as the input for our estimate. Another 
input value that favors the increase of the multiplicity is a relatively large 
value of the parameter a. This leads to larger values of the mean transverse 
momentum of emitted particles and of the integrated elliptic flow. 

We hope that our results will be useful for future comparisons between theory 
and experiments at the LHC. 
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